SUNY-NTG-96/1 



Universality near zero virtuality 



A.D. Jackson, M.K. §ener and J.J.M. Verbaarschot 

Department of Physics 
SUNY, Stony Brook, New York 11794 



Abstract 

In this paper we study a random matrix model with the chiral and flavor structure of 
the QCD Dirac operator and a temperature dependence given by the lowest Matsubara 
frequency. Using the supersymmetric method for random matrix theory, we obtain an 
exact, analytic expression for the average spectral density. In the large-n limit, the spectral 
density can be obtained from the solution to a cubic equation. This spectral density 
is non-zero in the vicinity of eigenvalue zero only for temperatures below the critical 
temperature of this model. Our main result is the demonstration that the microscopic 
limit of the spectral density is independent of temperature up to the critical temperature. 
This is due to a number of 'miraculous' cancellations. This result provides strong support 
for the conjecture that the microscopic spectral density is universal. In our derivation, we 
emphasize the symmetries of the partition function and show that this universal behavior 
is closely related to the existence of an invariant saddle-point manifold. 
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1 Introduction 



Since its initial application to level correlations in nuclear spectra QXJ] , random matrix the- 
ory has been applied to a variety of physical phenomena ranging from resonant cavities 
to lattice gauge theory [|J. In particular, the spectra of classically chaotic systems have 
been analyzed in great detail within this framework |4|, |5|, |6|. Recently, a great deal of 
progress has been made in the realm of mesoscopic physics. For example, universal con- 
ductance fluctuations have been understood thoroughly within the framework of random 
matrix theory 0, §, |, [TOf. 

In initial applications, random matrix theories were introduced in the hope of repre- 
senting the complicated strong interactions of nuclear physics. It was soon realized that 
the average spectral density cannot be described by means of random matrix theories. 
(For a large class of invariant random matrix theories, the spectral density is given by 
a semicircle. In contrast, the spectral density increases with energy for typical physical 
systems.) However, random matrix theories proved themselves capable of providing a re- 
markably accurate description of the correlations between eigenvalues on the scale of the 
average level spacing. Apparently, some of the properties of random matrices are univer- 
sal while others are not. Such behavior is familiar from the theory of critical phenomena 
where, for example, critical exponents are universal, but the critical temperature is not. 
There, universal phenomena are usually associated with the soft modes which arise due 
to the spontaneous breaking of a symmetry; non-universal properties are determined by 
all modes. A similar separation of scales takes place in random matrix theories. Consider, 
for example, the average spectral density (with variations only over many level spacings) 
and level correlations on the scale of the average level spacing. 

In the supersymmetric formulation of random matrix theory, universal properties, e.g., 
the level correlations, are associated with the existence of a saddle-point manifold which 
is intimately related to the symmetries of the theory. Non-universal properties, such as 
the average spectral density, can be calculated by a saddle-point approximation. 

In this paper we study a model which was introduced in [IT , O . This model is a 



random matrix model which possesses the chiral and flavor structure of the QCD Dirac 
operator and a schematic temperature dependence corresponding to the lowest Matsubara 
frequency. Otherwise, all matrix elements of the Dirac operator are completely random. 



The temperature dependence is such that this model has a second-order phase transition 
with mean field critical exponents. Below the critical temperature, chiral symmetry is 
broken spontaneously; above the critical temperature, it is restored. According to the 
Banks-Casher formula |[13J], the order parameter is the spectral density at eigenvalue zero. 



Physical motivation for this model comes from two rather different directions. First, at 
zero temperature, it satisfies all Leutwyler-Smilga sum-rules |15[], which are identities for 
chiral QCD in a finite volume. Second, Kocic and Kogut have recently suggested Jl6| that 



the chiral phase transition in fermionic systems is driven towards a mean field description 
because of the fact that the lowest Matsubara frequency non-zero. 

The chiral structure of the Dirac operator forces all eigenvalues to appear in pairs ±A. 
The spectrum is symmetric about zero. As we shall soon see, it is useful to introduce 
the microscopic limit of the spectral density which probes the spectrum around zero on a 
scale set by the distance between adjacent eigenvalues: 

1 u 

ps{u) = ^NE P{ NE ] ■ (1) 
Here, H is the temperature-dependent chiral condensate which, according to the Banks- 



Casher formula [13], is given by 

- TP(0) 



(2) 



N ' 

and N is the total number of eigenvalues. 

The zero-temperature version of this model has been studied extensively in the litera- 



ture |T7|, [18], [LJ], [2(J, |23|, |2l|| . It is known as the Laguerre ensemble or the chiral Gaussian 
Unitary Ensemble (chGUE). Two types of universal behaviour are known to exist in chiral 
random matrix theories: Spectral correlations in the bulk of the spectrum are universal; 
the microscopic limit of the spectral density, just introduced, is universal. It has been 



shown [22, 23 1 that the chiral structure of the random matrix ensemble does not affect 
eigenvalue correlations in the bulk of the spectrum. Such level correlations have been 
observed both experimentally and numerically in many systems ]T|, |2|, |], [| . Further, an- 
alytic arguments have been presented |5|, |(J in favor of the universality of correlations in 
the bulk of the spectrum of classically chaotic systems. 



In p5| we conjectured that the microscopic limit of the spectral density is universal 



as well. The first argument in support of this conjecture came from instanton liquid 



calculations [2f|, where we were able to generate ensembles large enough to permit the 



calculation of the spectral density in the microscopic limit. A slightly less direct argument 
came from lattice QCD calculations of the dependence of the chiral condensate on the 



valence quark mass [^7], Another hint came from the work of the MIT group ||29|| , who 
studied the chGUE using the supersymmetric method. They found that the microscopic 
limit of the spectral density is determined by a saddle-point manifold associated with the 
spontaneous breaking of a symmetry. The first convincing analytical arguments in favor 
of this conjecture came from recent work by Brezin, Hikami and Zee [pTf] . They considered 
families of random matrix models all possessing the chiral structure of the Dirac operator. 
They discovered the same microscopic limit in all the models they investigated. 

In this paper, we offer further evidence in support of the universality of the microscopic 
limit of the spectral density. Specifically, we investigate the effect of temperature in the 
chiral random matrix model introduced in O, 0|. This model differs structurally from 



the models in In the random matrix models considered in ||21|| , the unitary symme- 
try of the probability distribution leaves the spectrum of each element of the ensemble 
invariant. This invariance is not realized for temperatures T ^ 0, and analytic proofs 
are consqeuently somewhat more difficult. Using the supersymmetric method of random 
matrix theory, we obtain an exact expression for the spectral density which is valid for any 
dimension, n, of the matrices. This enables us to take the microscopic limit. This limit 
also requires the large- n limit of the spectral density (see (0)), which can be evaluated 
conveniently by means of a saddle-point approximation. In , this spectral density was 
evaluated numerically. It was found to have the well-known semi-circular shape at zero 
temperature. At high temperature, the shape is given by two disjoint semi-circles with 
centers located at ±7rT. In that paper, we also announced the analytic result for the 
shape of the average spectral density. The result merely requires the solution of a cubic 
equation. This result has also been obtained by Stephanov |5DJ , who applied an extension 
of this model to the problem of the relation between the Zn phase of the theory and the 
restoration of chiral symmetry. 

The organization of this paper is a follows. In section 2 we give a definition of the 
random matrix model and the supersymmetric partition function. In section 3 the par- 
tition function is reduced to a finite-dimensional integral. Symmetries and convergence 



questions of the partition function are analyzed in section 4. The exact two-dimensional 
integral for the resolvent is obtained in section 5. In section 6 we derive the large-n limit 
of the spectral density and discuss its properties. The microscopic limit of the partition 
function is evaluated in section 8, and concluding remarks are made in section 9. Our 
notation and conventions are explained in Appendix A, and a perturbative calculation of 
the large-n limit of the spectral density is presented in Appendix B. 



2 Definition of the random matrix model 



In this paper we study the spectrum of the ensemble of matrices 

( W + ttT] 
H -{w^ + nT ) ■ (3) 

Here, T is the temperature dependence as given by the lowest Matsubara frequency, and 
W is a complex n x n matrix distributed according to 

exv[-nY?TrWW ] ] . (4) 

The average spectral density can be expressed as 

2n 

p(A) = - lim — Im G(X + it) . (5) 



where the average resolvent G(z), 



2n z + iO-H 2n dJ 

can be obtained from the partition function 



(6) 



j=o 



Z(J) = [VW - d ^ T H \ T , exprflrW 1 ]. 
v ; J det(z + J-H) FL J 



(7) 



The integration measure, VW, is the Haar measure normalized so that Z(0) = I. For a 
Hermitean matrix, H, the resolvent is analytic in z in the upper complex half-plane. This 
allows us to calculate the resolvent for purely imaginary z and to perform the analytic 



continuation to real z at the end of the calculation. As will be seen below, this improves 
the convergence properties of the integrals in the partition function. 

Some properties of this model all already known. At T = 0, this model reduces to 
the well-known Laguerre ensemble. The joint probability distribution of the eigenvalues 
is known explicitly, as are all correlation functions. In particular, the average spectral 
density is a semicircle: 



We wish to stress that the largest eigenvalue is larger than a typical matrix element by a 
factor on the order of \fn. 



The temperature dependence of this model was analyzed in [11|. It was shown that, 
in the thermodynamic limit, this model shows a chiral phase transition at a critical tem- 
perature of 

Tc= ^ ■ (9) 

The order parameter is the chiral condensate H with H 7^ below T c and 2 = above 
T c . This chiral symmetry is broken spontaneously. For each finite value of n, H = 0. A 
non-zero value of 5 is obtained only in the thermodynamic limit. Below T c it was found 
that in this limit 

H = E (1 - tt 2 T 2 £ 2 ) 1/2 . (10) 

The Banks-Casher formula @ allows us to convert this value of S into the spectral density 
p(0). In [IT |, the complete spectral density of this model was determined numerically. At 



T = the result (H) was reproduced; at T = T c we found that p(X) ~ A 1 / 3 . For T ^> T c , 
the spectral density reduced to two semi-circles centered at ±7rT with a radius independent 
of T. We will present an analytic derivation of these results. 



6 



3 Ensemble average of the partition function 



In order to perform the Gaussian integrals, we write the determinant as an integral over 
the fermionic variables x an d X*'- 



U d [XiMXu]I[d[x2Mx2 



Similarly, the inverse determinant can be written as an integral over the bosonic variables 
d> and d>*: 



i „ n n 



The conventions for the Gaussian integrals are defined in Appendix A. The factor i in the 
exponent in (|T2|) is chosen so that the integral is convergent for z in the upper complex 
half-plane. This choice is consistent with the ie prescription in (|5p. The integral in ( |TT| ) 
converges independent of the overall phase of the exponent. The present choice of phase 
ensures that the product of the fermionic and bosonic integrals is one. 

The Gaussian integral over W can be performed by completing the squares. The result 
is a term of fourth order in the integration variables, 

exp --^{xljXu + 4>*2j(l>ii){XiiX2j + <j>*i<hj) • (13) 

We apply the Hubbard-Stratonovitch transformation to each of the terms of fourth order 
in the bosonic and fermionic variables. Two of the four factors can be decoupled with the 
help of real integration variables: 

exp \^4>\ • 0i 02 " 02 = / dai ^ a2 exp[-nS 2 (ai + o\) - {a x + ia 2 )0t • 0i + {a x - ia 2 )(p* 2 ■ 2 ] , 

7/.V ./ /;, 

1 * * f dp\(ip2 r 2 2\ / • \ * / 

exp +-=jXi -XiX2-X2= / — t — exp[-nE (p 1 + p 2 ) - (p 1 - ip 2 )Xi ■ Xi ~ (Pi + ip2)X2 • X2I ■ 

(14) 



7 



The terms which involve mixed bilinears can be decoupled with the help of Grassmann 
integrations. We do not encounter convergence problems in the process. 



( -i> -Z^,Xl • 01 X2 • 02 = 1 



exp 



nE 2 

1 



(Lot (Lj3 r o * r> * * ; n i * 1 

/ / "77/2)" exp ^~ nS " P + a Xi • 0i - P0 2 • X2J 



Xl • 01 X*2 ■ 02 



J/ / d (i/2) ex P[ -nS2 ^* a ~~ a *i ' 0i + ^*02 ■ X2] • ( 15 ) 



defined such that If/h = 1- (See Appendix A.] 



The constants If, and are 

Thus, we obtain the partition function as 

-nSVi + o\ + p\ + p 2 , + a*/? + /To) 
-ttT 



/n 
Y[ d[(f)i] d[xi\ d[a] exp 



x exp i 



/ 0* \ / Z + J + 20"i — (T 2 — 7T'J ' 2« 

02 — 7tT z + J — lO\ — cr 2 

X* i a * 2; + zpi + P2 

\X*2 J \ 2/3* -TTT 2 + Z>! - p 2 / V X2 / 





z/3 

-TlT 



\ /01\ 

2 

Xi 



(16) 



where 



d[cr] = fieri d&2 dpi dp2 



da da* d/3 d(3* 
WW 



(17) 



If ia\, o"2, z, and J are all real, the matrix A appearing in the exponent of (|T6|) is a graded 
Hermitean matrix. Then the Gaussian integrals can be performed according to (|73|) . This 
results in 



Z(J) = J d[a) exp -n£ 2 (a 2 + aj + pj + pj + a* (3 + /3*a) 



detg-M 



18) 



where detgA is the graded determinant of A. For a matrix with Grassmann blocks p and 
a and commuting blocks a and 6, it can be shown that 



detg ^ ° & ) = ^ ^ ^ e ^ a ~~ °^ 1/0 ^ 



(19) 



In our case, all blocks a and 6 are 2x2 matrices, which permits us to evaluate all 
expressions directly. 

We note, however, that the result (|TJ) was obtained by interchanging the 0j and <jj 
integrations in (|16|) . This is allowed only if the integral is uniformly convergent in a. 
Unfortunately, this is not the case when the 0\ and o<i integration paths are along the 
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real axis. This problem can be circumvented by a suitable deformation of the integration 
paths. Previous studies of random matrix theories within the framework of the sigma 



model formulation of the Anderson model |14j have stressed the importance of deforming 
the integration contours in a manner which is consistent with the symmetries of the 
problem. The same is true for the present problem. In order to interchange the and 



a integrals in (0), we must deform the integration contour so that the integration 
is uniformly convergent in a. In order to motivate our choice of contour, we must first 
consider the symmetries of the partition function. 



4 Symmetries 

We wish to study the partition function in the microscopic limit, i.e., the limit n — > oo 
with zn held constant. For z = (and J = 0), the partition function has an additional 
symmetry. For all temperatures, the bosonic part of the partition function is invariant 
under the non- compact symmetry operation 

01 - e +t 0i , 0t - e +t <P\ , 

2 - e-*0 2 , <Pl - e"V; • (20) 
This induces a hyperbolic rotation of the variables o~% and a 2 of the preceding section, 

(21) 



o~i \ I cosht isinht \ I °~i 
a 2 ) \ — i sinh t cosh t j \ a 2 



which clearly reveals the 0(1, 1) nature of the transformation. 
The fermionic part of the partition function is invariant under 



X2 - e-» X 2 , X*2 -> e- lu x* 2 , (22) 



where, a priori, u can be either real or complex. In terms of the p variables of fll4]) , this 
induces the transformation 



pi \ I cosn smti \ / p\ 
p 2 J I -sinw cosm J \ p 2 



(23) 



Because the integration over the Grassmann variables is finite, the volume of the symmetry 
group must be finite as well. Therefore, u must be real with u G [0, 2tt]. In other words, 
the symmetry group is 0(2). 

The terms of a mixed fermionic-bosonic nature are also affected by the transformations 



(p0|) and (|23|). This induces a transformation of Grassmann variables introduced through 
the Hubbard-Stratonovitch transformation. 

As is known from studies of random matrix theories, the parameterization of the vari- 
ables (<7i, (J2) and (pi, P2) is dictated by the above symmetries. It is natural and convenient 
to choose integration variables which lie along and perpendicular to the invariant mani- 
fold. In the microscopic limit, integrations along this manifold must be performed exactly, 
whereas the perpendicular integrations can be performed by saddle-point methods in the 
limit n — > 00. For the a variables, we thus choose the parametrization 

o\ = — i(a — ie) sinh s/E , 

a 2 = (z -ie) + J+ (a - ie) cosh s/E , (24) 

where a G [— 00; +00] and s G [—00; +00]. After the fa integration, the ie appears only 
in the combination o — ie. Below, we will not write the ie term explicitly, but it is always 
understood that it is included in the variable a. This parametrization renders the 0i and 
fa integrations uniformly convergent in o\ without jeopardizing the convergence of the a 
and s integrations. This allows us to interchange the fa and Oi integrations leading to the 
final result (|l^) of the last section. The term a\ + a\ appearing in the first exponent in 
(|i~8|) becomes a 2 + {z + J) 2 + 2(z + J)a cosh s in the parametrization (|24|). It is clear that 
the integral over s can be convergent only when z + J is purely imaginary. (Recall that 
cr contains the term — ie.) The transformation fl2"0|) for z + J = reduces to a translation 
of s and leaves a invariant. 

The p variables are also parametrized along and perpendicular to the saddle-point 
manifold according to 

Pi = iz + pcosip/T. , 

p 2 = psin<y9/E . (25) 

The rotation (|23| ) leads to a translation of the angle (p and leaves p invariant. 
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Our strategy in dealing with (jig) is to perform the Grassmann integrations first, i.e., to 
collect the coefficient of aa*(3f3*. This leaves us with a four-dimensional integral which is 
the exact analytical result for the partition function for any finite n. In the thermodynamic 
limit and for z ~ 0(1), the remaining integrations can be performed with a saddle-point 
approximation. This result is obtained in section 6. In the microscopic limit, z will be 
0(l/n), and the integration over the invariant manifold must be performed exactly. The 
radial integrals can be approximated to leading order in 1/n. (See section 5.) 

The T = problem has been investigated previously using the supersymmetric method 
[TjJ. In the saddle-point manifold was constructed in a manner similar to that 



used for the problem of invariant random matrix ensembles. (In this regard, see |3T], [32].) 
In |HJ] , the convergence difficulties were circumvented in an elegant fashion by the use of 
spherical coordinates for the variables </>i and <p2- Unfortunately, a direct generalization of 
this approach is not possible for the present case of non-zero temperatures. When T^O, 
the angles between the complex vectors <pi and 02 also enter in the integration variables. 



5 Exact result for the spectral density at finite n 



Because the fermionic blocks of the matrix A are nilpotent, the right side of (plj) can 
be expanded in a finite number of terms. The n-th power of the inverse of the graded 
determinant can be written as 



n 

a a 



dctg-"{" &) = (^) [l + ^Tra-VrV + ^Tr(a-VrV) 2 + yTr 2 a-VrV 



(26) 



Terms in the partition function which are of fourth order in the Grassmann variables 
can be obtained by supplementing the above terms with factors a* (3 and f3*a from the 
exponent in (|1|). The result is 



n 2 £ 4 f 
ix 1 J 



' tt 2 T 2 2 QD + 7r 2 T 2 )(A + 7r 2 r 2 ) +rr 2 T 2 (A- D)/n 
^ ~ DAT, 2 D 2 A 2 S 4 _ 

A\ n 

x (^J exp[-n£V 2 + a 2 + p 2 + p 2 )] , (27) 
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where D is the determinant of the boson-boson block, 



D = (z + J + i<J\ — cr 2 ) (z + J — io\ — cr 2 ) — 7T 2 T 2 



(28) 



and A is the determinant of the fermion-fermion block, 



A = (z + ipi + p 2 )(z + ipi - p 2 ) - 7r 2 T 2 



(29) 



In ( p7f) the variables cjj and p^ are parametrized according to ( ^4|) and fl25l) . tempera- 
ture by 



t = ttTS 



(30) 



the resulting form of Z(J) simplifies to 



Z(J) 



Ifl /"OO roc poo 

ada / ds 



7T- 



■2tt 

pdp I dip 
o 



f 2 



)2+ a 2 p 2 + t 2 (a 2 + p 2 )/n 



( t 2_ (T 2)( p2+t 2)/ ( (J 2_ t 2)2( p 2 +t 2)2 

-n(a 2 +p 2 +2(z+J)Sa cosh s+2jzEp cos ip+S 2 ((z+ J) 2 -z 2 )) 



> 2 +t 2 \ n 



(31) 



The integrations over s and <p can be expressed in terms of Bessel functions. Differentiation 
of the partition function with respect to J at J = gives us the resolvent which we desire. 
Thus, the final result of this section is 

2in 



G(z) = f odapdp 



X 



X 



t 2 )2 a 2 p 2 + t^ + p*)/n 

(t 2 -CT 2 )(p 2 +t 2 ) J (CT 2 -t 2 ) 2 (p 2 +t 2 ) 2 



y +t 2 ~ 

t 2 -a 2 



(ZzY^K^nllzo) - 2nZ(rK 1 (2ni:z(j))Jo(2ni:zp) exp[-n(cr 2 + p 2 )]. 

(32) 



Again, we remind the reader that a contains a term —is. The integral over a can thus 
be performed by successive partial integrations. Details regarding the calculation of this 
kind of integral can be found in ||33]| . This result has been obtained for z purely imaginary. 
Since the modified Bessel functions have a cut for za < 0, we can analytically continue 
this expression anywhere in the upper half-plane. 
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6 The large-n limit of the average spectral density 



For n — > oo and z ~ 0(1), all integrals in ( pT|) can be performed by a saddle-point ap- 
proximation. Because we started with a supersymmetric partition function, the Gaussian 
fluctuations about the saddle point give an overall constant of unity, i.e., Z(0) = 1 in 
(pl|). Using the relation (||) to determine the resolvent from the partition function fl3T|), 
we find that 

G(z) = £(£z + acoshs) , (33) 

where a and s are the saddle-point values of these variables. The saddle-point equation 
for s is trivial with solution s = 0. The equation for a is more complicated 

(a + Ez) = . (34) 



a 2 -t 2 

This equation can be rewritten as an equation for the ensemble averaged resolvent 

G 3 /S 4 - 2zG 2 /Y? + G(z 2 - tt 2 T 2 + 1 /S 2 ) -z = 0. (35) 
At T = 0, this equation reduces to 

(G - z)(G 2 /Y 2 - zG + I) = (36) 

with a non-trivial solution 

GW = E2 » ± i(4/g-^ (37) 

As indicated in (|) above, the associated spectral density is simply the imaginary part of 
the branch of G(z) with the negative sign, 

p (A) = !^!(4/S 2 -A 2 ) 1 / 2 , (38) 

7T 

which is the familiar semicircle normalized to the total number of eigenvalues. 
For z — 0, the saddle-point equation (|35|) simplifies to 

G 3 + S 4 G(1/S 2 - vr 2 T 2 ) = (39) 

with the solution 



G(0) = -iSVl - tt 2 T 2 S 2 . (40) 
13 



This corresponds to the spectral density 

On y 

p(0) = VI - 7T 2 T 2 S 2 . (41) 

TT 

Using the Banks-Casher formula (0), we immediately obtain the chiral condensate (|T0|) 



in agreement with [11 



In order to determine the high-temperature limit of the spectral density, it is most 
convenient to return to the saddle-point equation ([34]). It is clear that, for z ~ ttT, 
this equation can only be satisfied for a ~ —t. In the partition function (|3T| ) we can 
approximate the logarithmic term 

log(t 2 - a 2 ) « log(2t) + log(t + a) . (42) 

This leads us to the high-temperature limit of the saddle-point equation 

1 (a + Ez) = . (43) 



a + t 

The solution for the resolvent is 



G{z) = -pz-t-iy/2-pz-t)*) . (44) 

This results in a semicircular spectral density of radius a/2 located at z = txT. An identical 
argument leads to another semicircular contribution to the spectral density of radius a/2 
centered at z — —ttT. Of course, we can arrive at the same conclusion working directly 
from (|35|) . For z ~ ttT, the resolvent G(z) ~ 0(1), and the first term in (|35|) will be 
sub-leading in the high-temperature limit. This leads immediately to (|44). 



Finally, we consider the case at the critical temperature, T = E/7r, with z in the 
neighborhood of 0. Then, the saddle-point equation for G reduces to 

G 3 = z (45) 

with solutions (2S 4 ) 1 / 3 , (zE 4 ) 1 / 3 exp(7ri/3) and (2S 4 ) 1 / 3 exp(27rz/3). Only the last of these 
gives rise to a positive definite spectral density with 

p(A) = ^(AS)V 3 (46) 

7T 

in agreement with the mean field critical exponent of 5 = 3 for this model. 
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The equation for the resolvent (|35| ) also enables us to obtain a simple recursion relation 
for the moments of the spectral density. Expanding G(z) in terms of these moments, 



M- 



G{z) = E ~£l , (47) 

n 4/ 

we obtain 

M 2n+2 = (T 2 - ^)M 2n + — Y, M 2k M 2l -— Yl M 2k M 2l M 2m . (48) 

k+l=n k+l+m=n—l 

The evident initial condition, M = 1, immediately leads us to et cetera. Without too 
much effort, it is possible to use standard combinatoric methods to find the general result 
for the (2n)-th moment, 



7 The microscopic limit of the partition function 

The 'microscopic limit' denotes the investigation of the spectral density in the vicinity 
of z = on a scale set by the average level spacing. More precisely, we take the limit 
n — ► oo while keeping nz fixed, as indicated in (|1|). We start from the expression flT8"| ) 
for the partition function. In the thermodynamic limit, the a and p integrations can be 
performed by a saddle-point method. The saddle-point equations read 

p -p=o, 



p 2 + f< 

- a = , (50) 



t — u 
with solutions 

P 2 = l~t\ 

a 2 = t 2 -l. (51) 

For temperatures less than the critical temperature, p is real and a is purely imaginary. 
The integration range of p is the positive real axis. Therefore, the sign of p is positive. 
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The a integration ranges from — oo to +00. In order the reach the a saddle point, we 
must deform the integration contour. Because of the modified Bessel functions which 
appear in our expression (|32| ) for the resolvent, there is a cut in the complex a-plane for 
az on the negative real axis. The cut of the modified Bessel function is then ie above 
the positive real axis for negative z and ie above the negative real axis for positive z. 
Therefore, independent of the sign of z, only the saddle point with a negative imaginary 
part can be reached by a deformation of the contour. Thus, 

1/2 



a 



(52) 



for t < 1. 

At the saddle point, the pre-exponential factor vanishes: 

t 2 \ 2 a 2 p 2 + t 2 (a 2 + p 2 )/n 



1 



{t 2 - a 2 )(p 2 + t 2 ) 



{f -a 2 ) 2 (p 2 + t 2 ) 2 



. 



(53) 



Given (51), it is trivial that this equation is satisfied when t = 0. However, the vanishing 
of this pre-exponential factor for arbitrary t is remarkable and unexpected. This fact is 
responsible for the 'universal' behaviour of the microscopic limit of the spectral density. 
As a consequence, the 0(1/ n) term in this factor does not contribute to the resolvent 
to leading order in 1/n. The zKq term in the pre-exponent is also of subleading order 
(z ~ 0(l/n)). This leads to the following result for the resolvent in the microscopic limit: 



G(z) 



71 



X 



adapdp 
t 2 



+ 



a 2 p 2 



(t 2 -a 2 )(p 2 + t 2 )J (a 2 -t 2 ) 2 (p 2 + t 2 ) 2 
x aKi(2nEza) J a (2nLzp) exp[-n(a 2 + p 2 + log(t 2 - a 2 ) - log(p 2 + t 2 ))] . 



(54) 



In order to proceed, we make the substitution 

a = a + 5a 
P = p + Sp 



(55) 



in ( |54"|) and keep only those terms which contribute to leading order in 1/n, i.e., terms 
through second order in Sp and 5a. The exponent in (|54] ) then becomes 



exp(2na 2 6a 2 — 2np 2 5p 2 



(56) 
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The product of the terms in square brackets in ([54]) and a 2 p can be expanded as 

a 2 p[2p 3 5p + 2a 3 5a + p 2 (-l + 8t 2 )5p 2 + a 2 (+l + 8t 2 )5a 2 - ap{2 + 8t 2 )5p5a] . (57) 

It is clear already at this point that all temperature dependence enters through the scale 
factors a and p. The terms of 0(5p5a) vanish upon integration with the exponential 
factor (0). Since (Sp 2 ) = (5a 2 ) and a 2 = —p 2 , the other terms involving 8t 2 cancel as 
well. To complete the calculation, we need only expand the Bessel functions to first order 

K^nUza) J (2nJ:zp) = K x J + 2nEz8a K[ J + 2nHz5pK 1 J' Q . (58) 

where the Bessel functions K\ and Jo an d their derivatives appearing on the right of 
this equation are to be evaluated at their saddle points which are 2nLza and 2nY*zp, 
respectively. In order to arrive at the final result, we form the product of this expression 
and fl5TD, collect the coefficients of Sp 2 and 5a 2 , and perform the Gaussian integrations 
over Sp and 5a according to 

1 y/lT y/lT 



2 {2np 2 fl 2 (-2na 2 y/ 



2 



(V) 

[ ' 2 (2np 2 y/ 2 (-2na 2 ) 3 / 2 ' 1 ' 



The result is 



iT,a 2 r i 
G{z) = - — — KMa 2 - p 2 ) + Anzp^K^' + Anza^K^ J . (60) 
2 p A L J 

If we make use of the identities 

Jq = —Ji 1 

K[(z) = -K (z) --K x (z) , (61) 

we discover that the terms proportional to K\Jq cancel. This leaves us with 

G{z) = i2nzY?{\ - t 2 ){K x J x + iK J ) . (62) 

Finally, we can explicitly separate the resolvent into its real and imaginary parts by using 
two more elementary identities: 

M-iz) = —[h^+iN^z)} , 

K {-iz) = ^t[J (z)+tN (z)] . (63) 
17 



The final result for the microscopic spectral density is thus 

p(A) = 2n 2 A£ 2 (l - t 2 )(J 1 2 (2nASv / r 3 t 2 ") + J 2 (2nAS v / T 3 t 2 )). (64) 
As noted above, the microscopic limit of this model has previously been considered for 



the special case t = |25], |26], |29|, |2T|]. Our result is in agreement with this earlier work. 
Now, however, we can also consider the microscopic limit for general t 7^ 0. At finite tem- 
perature, the temperature enters only through the temperature-dependent modification 
of the chiral condensate which was obtained in |TT|. As defined in ([I]) with p(0) given by 
([|j]), the microscopic limit is strictly independent of the temperature: 

Ps {u) = -\Jl{u) + Jl{u)] . (65) 



8 Conclusions 

In this paper, we have studied a random matrix model with the chiral structure of the 
QCD Dirac operator and a temperature dependence characteristic of the lowest Matsub- 
ara frequency. This model possesses the global color and flavor symmetries of QCD. It 
undergoes a chiral phase transition with critical exponents given by mean field theory. 

Using the supersymmetric method for random matrix theories, we have found an exact, 
analytic expression for the average spectral density of this model. The result has the form 
of a two-dimensional integral which is valid for matrices of any dimension. In the large-n 
limit, these integrals can be performed using a saddle-point approximation. The spectral 
density then follows from the solution of an elementary cubic equation and nicely confirms 



our earlier numerical work 11 



Our primary result is that the spectral density in the microscopic limit is strictly 
independent of the temperature below the critical temperature of this model. This result 
supports the recent work of Brezin, Hikami and Zee, who investigated several families of 
random matrix models and found the same microscopic limit of the spectral density in all 
cases. As noted in the introduction, our model differs from the models considered by these 
authors in an essential way. In each of their models, the spectrum of each element in the 



ensemble is strictly invariant under the unitary symmetry of the probability distribution. 
In the present model this symmetry is violated for T 7^ 0. Thus, agreement between the 
microscopic limit of the spectral density for our model and the models of Brezin, Hikami 
and Zee increases our confidence in the universality of this quantity. 

In lattice QCD simulations the microscopic limit of the spectral density enters in 
the valence quark mass dependence of the chiral condensate. This quantity has been 



calculated for a variety of temperatures ||27|| , and it has been shown that the results below 
the critical temperature and not too large valence quark masses fall on a universal curve 
that can be obtained from the microscopic limit of the spectral density ||28|| . The present 
work provides a proper theoretical foundation of this analysis. 

In our derivations, the symmetries of the partition function played a crucial role. 
The universal behavior was closely related to the existence of an invariant saddle-point 
manifold generated by these symmetries. This suggests that the 'miraculous' cancellation 
of the temperature dependence of the microscopic spectral density found here is not 
a coincidence. It would be very interesting to obtain this result using more general 



arguments. Recent work by Guhr |3Jj on the superposition of two matrix ensembles 
appears to offer a promising method towards this goal. Work in this direction is in 
progress. 
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Appendix A: Notations and conventions 

In this appendix we summarize our notations and conventions. For a more detailed 
discussion regarding the motivation for these conventions, we refer to ]31, 32|. 



The integration measure for complex Gaussian integrals is defined such that 

exp +i<fi*<f) = 1 . (66) 



2vr 

For Grassmann integrals, the measure is defined so that 



2tt J dx* dxexp+ix*X = 1 • (67) 
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A graded vector or supervector is defined by 

-(J), (as) 

with <p a commuting vector of length n and x an anti-commuting vector of dimension m. 
The corresponding supermatrix which acts on this vector has the structure 

where a and b are complex matrices of dimension n x n and m x m, respectively. The 
entries in the n x m dimensional matrix a and the m x n dimensional matrix p are 
Grassmann variables. The graded trace of the matrix A is defined as 

TrgA = Tra - Tr6 . (70) 

The Hermitean conjugate of A is defined as 



where the f denotes transposition and complex conjugation. A graded matrix is called 
Hermitean if A^ = A. We use complex conjugation of the second kind for Grassmann 
variables, i.e., x** = ~X- The graded determinant is defined as 

detg A = exp(Trg log A) . (72) 

With this definition, we obtain the following natural result for a Hermitean, graded matrix: 

r n i 

J II d ^ rf M ex P = (73) 



Appendix B: Perturbative evaluation of the average spectral den- 
sity 

In this appendix, we derive the large-n limit of the resolvent without employing the 
super symmetric method. Because the operator (§) has only a finite support, it is possible 
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to expand the resolvent in a geometric series in \j{z — K) for z sufficiently large. Here, 
K is the matrix 

One finds by inspection that G(z) satisfies 



where is the matrix 



6 = 7^h- < 76 > 



and the bar denotes averaging over the probability distribution (f|). It should be clear 
that Q is block diagonal with the block structure 



gl n hl n 

hl n gin 



(77) 



where l n is the n x n identity matrix. Therefore, we find that G(z) = g. The average 
over W can be carried out immediately to give 



I W \ W \ 1 / gin \ 
^ H/t J y { J nS 2 ^ gin) ' 

This yields the following matrix equation for g and h: 

{ z -*T\(gh\l_(gO\(gh 
[-nT z )[h g ) + ^\0 g )[h g 

which leads to the two independent equations 

zg - TxTh = 1 + ^g 2 , 



(78) 



(79) 



zh-nTg = ±-gh. (80) 



S 



Elimination of h yields the equation 



irT*9 1 2 

zg ~ T^fs? = + t?' j ' < 81 > 

which agrees with fl34|) . Evidently, it can be rewritten as a cubic equation for g. 
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